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Universality of multicritical unitary matrix models is shown and a new scaling behavior 
is found in the microscopic region of the spectrum, which may be relevant for the low 



00 ' energy spectrum of the Dirac operator at the chiral phase transition. 
(N 

)Q ' 1. Introduction 
O 
00 

I In the last five years analytical results from random matrix theory have been successfully 
applied in the investigation of Dirac operator eigenvalues provided by QCD lattice data 
(see for a recent review). A crucial point for their applicability is the question of matrix 
Oh. model universality when replacing the Euclidean finite-volume QCD-partition function by 
r-| ! its matrix model counter part: 

Z = JdWdW^det'^f (^^ ^-NTrViWWf) ^ v{X') = j2^>^"' ■ (1) 

Here the chiral unitary ensemble with W N x [N + u) complex has been chosen as an 
example, which corresponds to QCD4 with Nf massless fermions, u zeromodes and gauge 
group SU {Nc > 3) in the fundamental representation. It is not fixed apriori which matrix 
model average exp(— A^V^) is to be taken since it cannot be derived from the effective 
Yang-Mills action, that averages the determinant of the Dirac operator in QCD. The 
matrix model calculation, if applicable, should therefore be independent of the details 
of the potential V and thus be universal. This question has been entirely answered for 
the unitary ensembles in the phase where the chiral symmetry is broken [0,0 . The link 
between the matrix model eigenvalue correlations and the chiral condensate is given by 
the Banks-Casher formula p(0) =| < qq >\/'n' [^]. However, the matrix model potential V 
can be also tuned in such a way that the eigenvalue density at the origin p(0) vanishes. 
These multicritical points may therefore serve as a class of models for the chiral phase 
transition, where the results presented here mainly summarize 0. In contrast to 
the transition is not driven by an external parameter introduced to model the effect of 
temperature T or finite baryon density /i in the matrix model. 



*Talk presented at the International Symposium QCD at Finite Baryon Density, April 27-30 1998 in 
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2. The chiral phase transition at multicriticahty 

There are different large- hmits in which the matrix model eq. (|lD can be investigated. 
In the macroscopic limit N{Xi — Xj) ^ 1, where the Aj are the eigenvalues of W, the 
eigenvalue density becomes a smooth function with finite support [—a, a] 

p(A) = J2Cki9^)X"'V^^^ . (2) 

fc=0 

It is nonuniversal since the coefficients Ckigi) depend explicitly on the couplings of the 
unknown potential V(A^) in the measure (see e.g. 0). Consequently the macroscopic 
density cannot lead to any quantitative prediction on the Dirac operator spectrum^. This 
has to be compared to the macroscopic density of the Dirac operator eigenvalues 

1 _ (iVi - 4) < gg >2 
pDiracW = <qq> + \^ F4 1^1 + ' 

which is clearly a non universal, model dependent quantity as well as it contains the pion 
decay constant physical parameter. 

Still, the macroscopic matrix model density carries qualitative information. If it is 
identified with the Dirac operator eigenvalue density |Tl| , the chiral condensate enters the 
matrix model through the Banks-Casher relation [Q. In the microscopic limit it will be 
the only physical parameter for all matrix model correlation functions which constitutes 
their predictive power [Q]. The fact that the first coefficient of the macroscopic density 
eq.(|^) C(){gi) ~< qq > is vanishing or not distinguishes between the matrix model being 
in the symmetric or broken phase. If the coupling constants Qi are adjusted such that the 
first m— 1 coefficients vanish, Co = ■■■ = Cm-i = (m-th multicriticahty), a class of possible 
models for the chiral phase transition is obtained, where the density vanishes as an even 
power p(A) ~ A^™ at the origin. A minimal potential describing such a multicritical 
behavior is given by 



kim) = 2^-^^^'^j^^-^'^-'^- , (4) 



(2m - 1) 

where the index + indicates to take only positive powers when expanding in 1/A^ and 
a = 1 has been chosen. The corresponding critical density reads 

pM) = ^A;(m)A2™v/r^ . (5) 
Ztc 

The parameter replacing < gg > is given by the first non vanishing term Cmidi) ~ p^'^™'\0) 
in eq.(0). The question of universality can then be addressed by perturbing the minimal 
potentials (|) by higher order terms while maintaining the same critical behavior. The 
phase transitions corresponding to these multicritical points are all of third order as their 
free energy F = l/A^^ln^ behaves like |Q 

F {g- 9.)'^^ , m e N , (6) 



^The connected density-density correlator (and all higher correlators) is universal in the macroscopic limit 
^ as it only depends on V through the support a. However, the support of the Dirac operator eigenvalues 
is a cut-off dependent quantity, such that macroscopic universality has no physical implications here. 
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when approaching the critical point. Investigations using matrix model P| or renormaliza- 
tion group techniques (see proceedings) have shown that the (T — /i)-phase diagram 
for < qq > consists of a first order (small T) and a second order line (small /i). In the 
tricritical point where the two lines meet the class of transitions provided by the above 
multicritical matrix models may become relevant. 



3. A new microscopic scaling limit 

In the microscopic limit correlations of eigenvalues are considered which are in a distance 
of the order Since the spectrum of the Dirac operator close to the origin is of interest 
a scaling variable x = NX is defined and kept finite as — oo. The rescaled microscopic 
eigenvalue density reads 

^ i^i?'^^^^^ = ^) • 
It is a universal function since it depends on the potential V only through the macroscopic 
density at zero, p(0), and thus it is parameterized by chiral condensate < qq > only. This 
has been shown for the unitary ensembles for an arbitrary polynomial potential at any 
given Nf and u for the massless P| and massive case P|. 

In order to extend universality to the multicritical points of the unitary ensembles the 
way to take the microscopic limit has to be modified. At the phase transition < qq >—>■ 
the appropriate scaling behavior for the m-th multicritical point is found to be 

N^^X = X , m e N , (8) 

and the microscopic density in terms of the new scaling variable reads 

Ps{x) = lim N~2^ p(^N-2^x) . (9) 

The same phenomenon may be expected for the Dirac operator eigenvalues on the lattice 
at the transition. The new universality classes eq.(|^) will be parameterized by p*-^™'*(0) 
in analogy to the broken phase 0. Recently the multicritical behavior of a different 



matrix model with random and deterministic degrees of freedom has been studied 
The authors find a modification of the microscopic scaling at multicriticality as well, with 

2fc+l 

exponents A 2^+2 A = x, A;gN. For k = l such a model has been applied to the chiral phase 
transition (see proceedings). 

4. Universality at multicritical points 

In this section a brief summary of the analytical and numerical results of is given. 
First the results for the unitary ensemble (QCD3) 

Z = /rfMdet2^/(M) e-^^^^(^) ~ T n Kl^^l'^^'e"^''^"^)) |det,,A}-f , (10) 

with potential V{M) = ^M'^^ are given and then extended to the chiral unitary en- 
semble. In order to determine correlation functions of eigenvalues a differential equation 
for the wavefunctions 

7/>^(A) = |A|^/e--^WP™(A) , / dX V^^(A)^„(A) = , (11) 
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is derived following [T^ and then solved in the microscopic hmit eq.(||). The P„(A) 
are polynomials orthonormal to the measure absorbed in the wavefunctions ipn{^)- The 
spectral kernel 

(I ,^ ^ V^jv(/x)^jv-i(A) - ^jv(A)V;jv-i(/i) 

H — A 

then determines all n-point correlation functions from taking its determinant. The set of 
polynomials P„(A) obeys the following properties 

AP„(A) = c„+iP„+i(A) + c„P„_i(A) , 

P„'(A) = A„(A)P„_i(A) - P„(A)P„(A) , (13) 

where the boundary condition determining the recursion coefficients c„ and the func- 
tions An{X) and Bn{X) can be found in |TB|. From eqs.(T3) the authors derive that the 



wavefunctions satisfy the following differential equation for finite 

C(A)-P;v(A)^;(A) + G^(A)^„(A) = 0, (14) 

where the PAr(A) and G]\f{X) are given functions of Aj^, Bjq and the potential V as well as 
derivatives of them |l^. Taking the scaling limit in the broken phase together with the 

relation p(A) =lim7v^oo ^'nv^ \J'^ ~ (-^Z*^)^ ([HI]) leads to a universal differential equation 
for of Bessel type first derived in When going to the m-th multicritical point 

in the scaling limit eq.(|^) the resulting differential equation will no longer be solvable 
analytically, the P/v(A) and Gj\f{X) being rational functions 0]. 

In the simplest example the m = 1-critical potential V{X) = — 4A^ + ^gX"^ at (7 = 16 leads 
to the following coefficients in terms of the scaling variable x = N^^^X: 

iV-lp^(x) = (15) 

The constants u± = g{2p{{]) ± /'(O)) and v = 2gf{0) have to be determined numerically 
by solving an auxiliary Painleve II equation for the recursion coefficients at ^ = 



= gf{z)' - zm - p"{z) + ^ . (16) 

Still, the issue of universality of the two differential equations and thus of ipNi^) can be 
addressed analytically. Perturbing the critical quartic potential by a sextic term geX^/G 
and maintaining p(A) ~ A^ all expressions eqs.(|T5|) and (|1^) remain valid when replacing 
g ^ g* = g + g&l2. Thus g^, = p"{0)/n plays the role of a universal parameter in the 
solution iPn{x) and via the kernel eq.(|12|) in all microscopic correlation functions. Below 
the numerical solution for the universal microscopic density eq.(^ for m = 1 is given. The 
x^-growth is due to the matching condition lim^^^ pg{x) = p" (0) (dotted line). 
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Figure 1. The m = l multicritical density at Nf = (left) and Nf = l (right). 

For higher m the functions F]sf{x) and Gn{x) become much more involved. However, 
when making an approximation to the exact differential equation (III) in the scaling limit 
(P) it is possible to obtain an equation which is analytically solvable for any m. Imposing 

1 < k{m) x'™ (17) 

on X for the critical potentials (^) leads to the approximate equation 

X \2[2mlj x'^ J 

which is again of Bessel type. Here the constant A;(m) = 27rp*^^'"^(0)/(2m)! has been 
replaced. The solution read^ 

J2m)(n\ 

^^{x) ~ VXJJ^ ^{X) , X= P \l x"-^' . (19) 
2m+i 2 \zxn -r ij! 

The appearance of the parameter p'^^™^(0) in the approximate solutions makes it highly 
suggestive to conjecture new universality classes for all multicritical points with m>2 . 

The above results can be easily translated to the multicritical chiral unitary ensembles. 
The wavefunctions '?/'n(A) corresponding to the partition function eq.(|l|) can be simply 
related to those of the unitary ensemble above 

^^{\^- iVj + Z/) = i)2m{\ iV^ + Z/+^) , (20) 

when shifting iVj in eq.(p!0|) by z/ + 1/2. The approximate analytic solution can thus be 
immediately read off from eq . (p!9D . The final expression for the approximate microscopic 
density of the m-th multicritical chiral unitary ensemble reads 



_ „(2m) /rj\„2m / 9 A/^ 



with P = -^^^(2m^^- matches to the exact noncritical density for m = 0. 

•^To fix a unique solution regularity and normalizability have been imposed as in the broken phase, 
although eq.(18) is no longer valid at x = 0. 
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Figure 2. A comparison between the numerical solution for the m = 1 critical micro- 
scopic density for Nf + = (upper line) and the analytical approximation (lower line) 
which breaks down at small x (see eg. (P^))- 
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